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Abstract 

We study in detail the transport properties of a model of conducting electrons in the presence of 
double-exchange between localized spins arranged on a 2D Kagome lattice, as introduced by Ohgushi, 
Murakami and Nagaosa (2000). The relationship between the canting angle of the spin texture 6 and 
the Berry phase field flux per triangular plaquette 4> is derived explicitly and we emphasize the 
similarities between this model and Haldane's honeycomb lattice version of the quantum Hall effect 
(Haldane, 1988). The quantization of the transverse (Hall) conductivity Oxy is derived explicitly from 
the Kubo formula and a direct calculation of the longitudinal conductivity a^x shows the existence 
of a metal- insulator transition as a function of the canting angle 9 (or flux density (f>). This transition 
might be linked to that observable in the manganite compounds or in the pyrochlore ones, as the spin 
ordering changes from ferromagnetic to canted. 

1 Introduction 

A fascinating and relatively common problem in solid-state physics is the motion of electrons in a 
lattice structure where localized spins are present, possibly in an organised form. The scattering of the 
electrons by the spins represents a complex physics problem for the theoretical investigator, especially 
when (like it seems to happen for the high-Tc cuprate superconductors) the spin and charge degrees 
of freedom are attached to the same particles - a situation which we do not consider. In a simplified 
picture we assume first of all that the spins are localised and interacting, but that the mobile electrons 
do not interact with each other (or interact weakly and give rise to independent quasiparticles) and are 
separate entities from the electrons that produce the magnetic ions. In this case the individual electrons 
(or quasiparticles) experience the localized spins as an effective local magnetic flux (Miiller-Hartmann 
and Dagotto, 1996). This flux produces a Berry (or Peierls) phase in the hopping terms of the electrons' 
Hamiltonian. In this paper we shall study a situation with a locally staggered flux, created by plaquettes of 
ordered spins, where the global Rux is zero. This case is realized, for instance, in a Kagome lattice (Ohgushi, 
Murakami and Nagaosa, 2000), a two dimensional (2D) lattice consisting of triangles and hexagons of 
the same interatomic distance and that can be viewed as a triangular Bravais lattice with a three-point 
atomic basis forming an equilateral triangle of size half that of the triangular lattice parameter's. In Fig.l 
we remind the reader of the Kagome lattice structure. 

This type of lattice may have experimental relevance in the planes of pyrochlore compounds (Ohgushi, 
Murakami and Nagaosa, 2000, Ramirez, 1994, Harris and Zinkin, 1996) and the transport properties we 
describe may be appropriate for such materials. Ferromagnetic pyrochlore crystals of the type R2M02O7 
(R=Nd, Sm or Gd) have revealed interesting transport properties like an anomalous Hall effect increasing 
as the temperature T is lowered, (Taguchi and Tokura, 1999) a feature that seems to be connected with the 
geometrical frustration of pyrochlore lattices that is partly embodied by the Kagome lattice itself, viewed 
now as the (1,1,1) cross-section of the 3D pyrochlore's lattice. Motivated by the transport properties of 
pyrochlore compounds, as well as by those of the manganite ones, we study in this paper some quantum 
transport properties of the Kagome lattice with a canted localised spin texture in which independent 
electrons can move. 
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Figure 1: Kagome lattice 

The behavior of the Berry phase as a function of the spin canting for the model at hand, and its 
consequences in terms of the macroscopic transport properties of the model are studied in this paper. 
We discuss in detail the energy spectrum as it depends strongly on the canting of the localized spins. In 
particular, the nodes in the spectrum and the opening of energy gaps are investigated, including their 
consequences for the transport properties. Moreover, we evaluate explicitely the longitudinal conductivity 
(Txx and the Hall conductivity (Txy as a function of the canting angle 9 or flux cf). Our results for the 
transport properties are then compared with those of another famous model of this class, where a staggered 
magnetic field is applied to electrons within a honeycomb lattice (Haldane, 1988). The latter has very 
similar spectral properties as the model, first proposed by Ohgushi, Murakami and Nagaosa, (2000) 
defined on the Kagome lattice and studied here in greater detail. The model on the honeycomb lattice 
was proposed by Haldane as the condensed-matter (solid-state) equivalent of the quantum Hall effect, in 
that a quantization of the Hall (Txy conductivity can be achieved hy varying the local flux per plaquette 
(/>, but without the need to introduce an external, homogeneous magnetic field. In the case of the Kagome 
lattice model at hand, the same result will be shown to be attained through the introduction of a localised 
spin texture. Our results confirm and complete the work by Ohgushi, Murakami and Nagaosa for the 
Hall ffxy conductivity, with exphcit calculations in terms of expansions around the gap's nodes shown in 
detail, and moreover the quantized values of the longitudinal Uxx conductivity are obtained indicating the 
existence of some sort of metal-insulator transition as the canting angle moves away from some special 
values. 

The paper is organized as follows: In Section 2 the tight-binding model of Ohgushi, Murakami and 
Nagaosa for localized spins is described. The relationship between the localized spin's wavefunction and 
the Berry phase of the hopping elctron is discussed in Section 2.1, and the apphcation to the Kagome 
lattice (Section 2.2) is presented. In Section 2.3 the tigh-binding model for the honeycomb lattice with a 
staggered magnetic field is also presented and compared to the model on the Kagome lattice. Transport 
properties are then studied in Section 3 by making explicit use of Kubo's formula and an expansion of 
the energy spectrum near the nodes next to the Fermi energy, and the results obtained are discussed in 
Section 4. 



2 



2 A Model of Hopping Electrons in a Spin Texture 
2.1 The Model 

We consider the electronic hopping between nearest neighbours on a Kagome lattice as described 
by the tight-binding Hamiltonian. The electronic degrees of freedom are coupled to a set of localized 
spin-S degrees of freedom on the same lattice via a local Hund coupHng Jh- In this work S = \, but 
generaHzation to the physically and theoretically interesting case of larger S is possible. When Jh is 
strong enough the spin of the hopping electron is forced to allign parallel to the localised spin Sj at each 
site and through double-exchange mechanism (Zcncr, 1951, Anderson and Hascgawa, 1955, dc Gcnncs, 
1960 ) the tight-binding hopping parameter Uj between two neighbouring sites becomes proportional 
to the projection of the localised-spin wave function at site j onto that at site i. The effective Hamiltonian 
representing the hopping is then 

H=Y.ni^cic,+h.c. (1) 

where tl^^ = t(nj|nj), t being the bare hopping parameter and |n) the spin wave function for a spin-^ 
quantized along the direction defined by the unit vector n = (sin 9 cos 0, sin 9 sin (j), cos 9) . This spinor 
wave function clearly satisfies (with a = {ox,Oy, Uz) the vector of Pauli matrices) n • CT|n) = -|-|n) and is 
given by 

I") = ^"(e52lf )' 

where h is an undetermined overall gauge degree of freedom. The effective hopping parameter is then 

tiff = te-'^b,-'^^ (cos 5 cos ^ + e-^(^*-^^) sin | sin ^| (3) 
[^2 2 2 2 J 

and since |(nj|nj)p = cos^ with cos 6ij = • nj or Oij being the angle between the two localized 
spins' directions of quantization so that cos^ ^ = ^(1 -|- iij • nj), we see that we can put 

tiff = fe»« cos ^ (4) 
where the Berry phase Uij is obtained (ignoring the gauge parameters) through 

cos % cos % -I- e-'('^*-'^^) sin sin % 

giay ^ 2 2 2 2_ (^g-j 



sin Y sin -f- sin(0j — (pj) 



and can be evaluated, e.g., by means of 

sin aij = - T'"^ • (6) 
cos -f- 

To see what the phase Uij is, geometrically, we introduce the unit vector z and evaluate the triple product 

9' 9 - 9- 9 - 

Ui X Uj ■ z = sinOi sinOj sin((^i — <pj) = 4sin sin cos cos sm{<pi — (7) 
which shows that 

X n,- • z 

^^"^ = ^ J e- ■ (8) 

4 cos -f cos -J cos -if 

This expression is a special case of the formula giving the solid angle f2(ni,n2,n3) between three unit 
vectors ni, n2 and ns: 

. 0(ni,n2,n3) ni x n2 • na 

sin 7, = 3 3 5— (9) 

2 4 cos % cos % cos % 
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(as can be verified by taking, e.g., ni — x, 112 — y and ns = z; fl{ni, n2, n^) can also be seen as the area 
of the portion of unit sphere enclosed by the maximum circles passing through the unit vectors' tips). In 
the last formula, of course, cos 9kk' = nfc • n^/ and therefore cos^ = ^(1 + • n^/). We remark that 
this formula for three spins is completely analogous to that for the chirality gauge field in the formulation 
of Lee and Nagaosa (1992) for the chiral spin hquid theory of high-temperature superconductivity (Wen, 
Wilczek and Zee, 1989). In this formulation the instantaneous gauge field fiux through the triangular 
plaquette made up by the three spins is <i>(ni, 02, na) = iri(ni, n2, na). Back to our two-spins hopping 
problem, we then conclude that sinoij = — sin iri(ni, rij, z), or 

a^j = TT +^n{ni,nj,z) (10) 

(the factor ^ being probably due to our specially chosen localised spin value, which leads to the conjecture 
that for a generic spin-S* situation the Berry phase would be = tt -I- 5^2(11^,11^,2)). Since the solid 
angle fl{ni,nj, z) is also the unit sphere's surface area between the tips of the three vectors n^, rij and z, 
the phase = dr ■ A can also be seen as the fiux of a magnetic monopole's field of modulus |B| = ^ 
with the monopole placed in the sphere's center, or, alternatively, as the fiux of the related gauge field 
A through the triangle bearing on the segment of a triangular Kagome lattice's unit cell (Ohgushi, 
Murakami and Nagaosa, 2000). In this way, the Berry phase acquires some physical meaning too. 

We now consider this tight-binding model on the Kagome lattice with a fixed localised-spin config- 
uration (or spin texture) as suggested by Ohgushi, Murakami and Nagaosa (2000) , in which the unit 
vectors at each site of a triangular unit cell are tilted outwards at a fixed angle over the unit vector 
z orthogonal to the lattice plane. This means (labelling the spins clockwise A, B and C in the unit cell) 

sma^B = smasc = smacA = . (11) 



4(1 + cos6i)y'l - ^sin^e 
The fiux generated by the spins in every triangular unit is set equal to ip with the condition 

gi0 _ gi{a.AB+a,BC+acA) _ gSia^s (12) 

with (j) = Sqab (mod 2tt) and thus 

. V3(l-cosg) 

sm- = — =. (13) 

3 2Vl + 3cos2 6i ^ ^ 

This is equivalent to the expression proposed by Ohgushi, Murakami and Nagaosa (preprint of Ohgushi, 
Murakami and Nagaosa (2000)) <f> — tt + 3arg (1 - iVS cos 9). The graph for this expression of = 0(0) 
is shown in Fig. 2 for convenience. 

As pointed out by Ohgushi, Murakami and Nagaosa (2000), the fiux per triangular unit cell (j> is 
cancelled out for the chosen spin texture by the fiux —2(/) generated by each of the remaining hexagonal 
hopping plaquettes on the Kagome lattice. There are indeed twice as many triangular units as hexagonal 
plaquettes, so that the overall gauge field fiux through the lattice is zero. This situation is reminiscent 
of the analogous tight-binding model in a staggered magnetic field as was proposed by Haldane (1988) 
to mimick the quantized Hall effect in a condensed-matter situation. In the present model, the chosen 
spin-texture, with all localised spins tilted by the same angle 9 in each unit cell, is presumably the 
one corresponding to the mean-field solution for some magnetic spin-spin Heisenberg Hamiltonian which 
should be added to our tight-binding Hamiltonian, Eq. (QJ, to give a total Hamiltonian of the type 

Htot = E i^/^({sf ^})clc, + h.c. + E -/..S, . S,. (14) 

The role of the spin fiuctuations around this ordered spin texture, {Sj-*'''}, as well as the effects of different 
spin textures, (e.g. AFM ones) could serve as an interesting further research problem for future studies. 
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Figure 2: Flux </) as a function of the tilting angle 9 of the localized spins. 

2.2 Band structure for the Kagome lattice 

The Kagome lattice is made up of triangular and hexagonal plaquettes and can be seen as a triangular 
lattice with a 3-point basis where every triangular unit cell contains three sites, A, B, C (Fig 1). The 
displacement vectors between these sites are ai = (— 1/2, — V3/2), 0.2 — (1,0) and 03 = (—1/2, a/3/2), 
with J2i — 0- The effective hopping parameter of 1^ can be written as 

i-f =<e»-cos(^) (15) 



and since cos(&y /2) = y 1 — | sin 9 is fixed for the chosen spin texture, 6ij being the angle between 

the n.n. pair of localized spins, we can choose the convention where tcos(^) = 1. Then, in momentum 
space, the Hamiltonian can be rewritten as 

H ^^ijj\k)h(k)ijj{k) (16) 
k 

where ij^ik) = (^CA(k), csik), cc(k)^ and h{k) is the (suitably symmetrized) matrix 

2cos(A?- ai)e^*'^/3 2 cos(fc • a3)e*'^/3 

h{k) ^ 2cos(fc-ai)e*'^/3 2 cos(fc • 02)6-*"^/^ |. (17) 

2cos(fc • a3)e-*'^/3 2 cos(fc • a2)e*"^/3 

The three eigenvalues of this Hamiltonian are 



-B™d(fc) = 41/ ^cos(^^- ) (18) 



,T^ . + ,0(A:) + 27r, 
Edownik) = 41/ ^^cos(^^^- ) 
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with 



0{k) = arg 



f{k) cos((/)) + i 



\ 



1 + f{k) 



(/(fc)cos(0))2 



(19) 



and 



f{k) = 2 cos(A; • ai) cos(A; • 02) cos(fc • 03). 

The three bands touch in six points only for (6 = and (f) = ±ir, while for 
there is a gap between the bands, as shown in Fig. 3 = 0) and Fig. 4 (</> ■ 



(20) 

different from these values 
:7r/3). 



E(k) 




Figure 3: Eigenvalues for the Kagome lattice for ^ = 

E(k) 




Figure 4: Eigenvalues for the Kagome lattice for (p — 7r/3 

These nodes are: (±^,0) and (±^, ±^7r), on the vertices of a hexagon, as shown in Fig. 5. 

The problem of calculating transport coefficients with this 3x3 matrix is not exactly solvable, so 
we reduce this matrix to a 2 x 2 one by expanding h{k) around the nodes. This can be done with a 
unitary transformation which allows us to neglect the terms related to the lower band; in fact this band 
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Figure 5: Nodes position for the Kagome lattice 



is far from the other two and gives no relevant contribution to the Green function present in the Kubo 
formula. To find this unitary transformation we consider the Hamiltonian evaluated at a node kp. If we 
apply the unitary matrix that diagonalizes h{kxo, kyo), where {kxo, kyo) are the node's coordinates, to the 
Hamiltonian evaluated at the general point {kx,ky), we find a matrix H' with the structure 



H' = [ a[kx+ iky] 
P[ka; + iky] 



a[kx - iky] 
7[fcx - iky] 



P[kx iky] 
•jlkx + iky] 
As 



(21) 



with kx = {kx — kxo) and ky = {ky — kyo). The elements on the diagonal are the eigenvalues of the 
Hamiltonian h{kxO;kyo), while the ofF-diagonal elements are complex combinations of fc^.y. Near the 
nodes the distance between the upper and the middle band is small, while the lower band is distant and 
gives no relevant contribution. To justify this we can consider a projection of the Green function around 
the node. 

We want to find a projection only on the first two eigenvalues, so we choose a projector such that 



PH P = 



Ai 

H21 ^2 



(22) 



We define the Green function as G = {z — H and the projection operator P with the convention that 
{A)p^ = {PAP)~^ is the inverse operation on the projected space, (1 — P) being the projection operator 
complementary to P. The Green function can now be written as 

G = PGP + (1 - P)GP + PG{1 - P) + (1 - P)G(1 - P) (23) 

and the projected Green function is 

PGP = P{z - h')-^P ={z- PH'p - PH'{1 - P){z - H')-^p{l - P)H' P)'^ . 

The last terms are of higher order and can be neglected. If we consider the terms related to the lower 
eigenvalue, we can see that 1/(2; — A3) ^ 1/2, because we are considering k near the nodes. This eigenvalue 
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does not give an important contribution, so it can be neglected and we can use the approximation 

PGP K {z - PH' P)-^ . (24) 
Now we can write the projection of the Hamiltonian H' as 

PH'P = ( ~ '^'^ ~ "[^^ ~ ^^v\ 
a\kx + iky] A2 



(25) 



/ 1 + hi^ih2\ 

\ hi+ih2 1-4^ y 

\ ni+ in2 —m j 

where I is the identity matrix, with 

hi = aikx + oi2ky 

h2 = -a2kx + aiky (26) 
m = —(j), 

with ai and 0:2 are the components of a complex parameter, depending on the node we are considering. 
The new matrix representing the Hamiltonian has eigenvalues ±A and eigenvectors ^± with 



A = -^TO^ + hl+hl 
1 



1 



Vl + (±A-m)V|fc|2 V ^ 

As was said, the bands touch only when </> = and (j) = iTi". Here we consider </> different from zero, but 
small enough so that we create a small gap between the two bands. 

2.3 Band structure of the honeycomb lattice 

Here we also consider the case of a honeycomb lattice, as was first envisaged by Haldane (1988); this 
is made up by two sublattices that we call A and B (Fig. 6), or by a triangular lattice with a 2-point 
basis. The displacement vectors from a B site to the three nearest neighbours are: ai = (— v/3/2, 1/2), 
02 = (0,-1) and 03 = (•\/3/2, 1/2), while the displacement vectors from the site B and the nearest 
neighbours on the same sublattice are bi = (-\/3/2,-3/2), 62 = (V3,0) and 63 = (-^3/2, 3/2) (again, 
E^.=0). 

Here too, we consider a tight-binding model in the presence of a staggered magnetic flux (Haldane, 
1988); the Hamiltonian for this Haldane model is H = {k)h{k)ilj{k) with 

h{k) = 2*2 cos (j) cos(k • bi)I + ti ^^[cos(k • ai)cri + sin(k • ai)(T2] 

i i 

+ [M-2i2sin(^^sin(k-bi)](T3, (27) 

i 

where ti is a hopping parameter between nearest neighbours on different sublattices, t2 is a hopping 
parameter between nearest neighbour sites on the same sublattice, and crj are the three Pauli matrices. 
If we rewrite the Hamiltonian as H = aai + h(j2 + ca^ , the Hamiltonian matrix reads 
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Figure 6: Honeycomb lattice 



The eigenvalues of this matrix are ±A with A = vc?+P + c? (Fig- 7 = 0) and Fig. 8 ((^ = ^)), 
while the eigenvectors are 

^± = ^=^f±L) (29) 



where k — a ~ ib. Formally this case is similar to that of the Kagome lattice, after the reduction of the 
original matrix to a 2 x 2 one. The two bands meet when the condition + 6^ + = is satisfied. This 
becomes a condition on the parameter M: there are nodes when M = 3-\/3i2asin0, with a = ±1. When 
= and Af = there are six nodes: (±^, 0) and (±^, (Fig. 9), while when M = 373*2 a sin </) 

there are only three of these nodes. 



3 Transport Properties 

Based on the Hnear-response theory, a suitable Kubo formula and the corresponding conductivity 
tensor can be studied for the Hamiltonians considered in Sect. 2. Some details are given in Appendix A. 
From this result we can derive the Hall conductivity axy and the longitudinal conductivity axx of our 
two-dimensional tight-binding model. We verify explicitely the quantization of function of (f) (or 

0) and calculate explicitely the longitudinal conductivity axx which also appears to be quantized in the 
absence of disorder or other symmetry-breaking conditions. For the two models considered, these are our 
main new results. 

3.1 Hall conductivity 

For ^ ^ V the third term in Eq. ll62)l (Appendix A) vanishes and after the integration with respect 
to E we find that the Hall conductivity, in the limit of w = and T = 0, is 

(T.y = -^i?e^ Afc [U^Jt,{h-\k + 2ir])-^j,U]kk, (30) 

^ k 

where U is the unitary matrix that diagonalizes the Hamiltonian matrix h{k), while jx and jy are the 
current matrices. From this we find that the Hall conductivity for every node n for the case of the Kagome 
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Figure 7: Eigenvalues for the honeycomb lattice for </> = 

E(k) 




Figure 8: Eigenvalues for the honeycomb lattice for (j)= ^, M ^ and M ^ 3\/3at2 sine 



lattice is 



"^^^ hr] 



+00 



12mr] (Pk _ 



(31) 



The integration boing ovor tho hoxagonal Brillouin zone, wc considered only one third of the integral and 
then we have to multiply for the number of nodes. We can conclude that the Hall conductivity is different 
from zero and it is quantized in the presence of a gap between the bands (that is for (f) different from 0, 
±7r) and is equal to 



<^xy = —sgn{(j)). 



(32) 



So, we have another model of transverse conductivity quantization in the absence of an external uniform 
magnetic field. 

Now we consider the case of the honeycomb lattice, in Haldane's model. For M = and = the 
bands touch in six points: (±^^,0), (±^^,±^) and the Hamiltonian is simply of the form 



a — ib 
a + ib 



(33) 
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Figure 9: Nodes position for the honeycomb lattice for = 

Expanding the terms around the nodes, we find that the function to integrate in order to find the Hall 
conductivity is 

± 3 (34) 

8(o2+&2)f 

but the integral of this term gives zero contribution. To generate a gap we have to move from the situation 
in which M = and <!) = 0. We add a small mass contribution M <^ 1, but we mantain = 0; now the 
Hamiltonian is 



M a-ib 
a + ib -M 



(35) 



In this case the function to integrate is 



± nMl^, (36) 

8(a2 + 62+M2)f ' 



three nodes give a positive contribution and three a negative one, to give 

± ^sgn{M) (37) 

Summing up all the contributions we find that the Hall conductivity is zero. Different is the situation in 
which we consider M = 0, but we add a small fiux 0. We rewrite the Hamiltonian as 



c a — ib 
a + ib —c 



(38) 



where 

c = -2t2 sin ^ sin(k • bi). (39) 
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Now the function to integrate is (fc = a — ib) 

P{k) = ± ^£!^!L^. (40) 

8(a2 + 62+^2)1 

Near three of the six nodes P{k) is negative, but approximating c around these points we find c ~ 
sin(/), so the function to integrate is negative. After the integration we find that every one of these 
three points gives a conductivity equal to 

—sgn{sm(j)). (41) 

For the other three points P{k) is positive, but the expansion of c is c ~ — 3\/3i2 sin(0). So, now the 
function to integrate is negative too and it gives the same result as before. Summing over all the points 
and remembering that the integration is over the exagon we find that the Hall conductivity for the 
honeycomb lattice case is 

CTi-y = — sgn(sin(/)). (42) 
We can conclude that the Hall conductivity can be rewritten as 

h 



crxv = v— (43) 



with V ~ ±1, depending on the sign of (j). 

3.2 Longitudinal conductivity 

Here too we use the general expression derived from the Kubo formula. In this case {fi = v = x) the 
longitudinal conductivity axx derived by expression Ij62|l . after the energy integration, is 



h ^-^ —uj + 2iri 



\k — uj + 2irj A„i — Afc + a; — 2irj 



After summing over the eigenvalues, substituting the values of Afe and e = 77 + ^ we find 



(44) 



>Prom this expression, we have still to subtract the diamagnetic term and so we evaluate (Ludwig, Fisher, 
Shankar and Grinstein, 1994) 



XX -limeaxx- (45) 



For the Kagome lattice the product of the matrix elements of the currents is 
(TT^ ■ n\ (TT^ ■ jn 3(4m2 + 3/1^ + 2^/3fel/l2 + fej) 

[UhxUUiUhxUh, = 4(^2 + . (46) 

and remembering that hi and /i2 are symmetric variables and using polar coordinates we can rewrite a^x 
as 



4he 



-dr + / — -dr 



%/to2 + r'^{m? + r2 + e^) Jq m? + r"^ + e 



(47) 



so that carrying out the integrals we find 



„ I ^ arccosf , 7, ■> ) if m > 

-dr= I Q if m = (48) 



i\rn? + r^{m? + + e^) 



-(tt — arccos( /"?,, ) if 



m < 
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For the second integral we introduce a cut-off A and evaluate 



r\J m? + r' 

-dr = < 



e(m-^ + + e^) 



|(-TO + Vm2 + A2) + arctan(2t) - arctan(^^^S2!) if m > 
|-arctan(4) if to = (49) 

^(m + + A2) - arctan(2i) - arctan(^^^SAl) if m < 0. 



We consider two cases: m equal to zero (that is, the flux (j) is zero) and m different from zero. In the flrst 
case, after having removed the diamagnetic term, the conductivity is 

e2 A. 

Cj;!; = ^ arctan(-). (50) 

Now we can take the limit for e ^ (thus making the cutoff irrelevant). We find that the conductivity 
for every node is different from zero and is equal to 

After the sum over all six nodes is done we can conclude that the longitudinal conductivity for to = is 

(52) 
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e TT 
'Ah' 



The case where to is not zero is quite different. Now, after having done the diamagnetic subtraction Ij45|l . 
the conductivity is 



e 



2 to2 to , to to V to2 + A2 



„^-^ arccos( , ) (- arctan( — ) — arctan( )1, (53) 

■ih' e Vm2 + e2 e ^e' ^ e " 

but the Hmit for e ^ gives a vanishing result. We conclude that the longitudinal conductivity is different 
from zero only when there is no gap between the two bands, that is in our case for to = 0. 

Next we consider the longitudinal conductivity for the honeycomb lattice. We expect that it is 
different from zero when the bands touch. This happens in six points, when M = and = 0, and 
in three points when M = 3\/3ai2 sin(0). To calculate axx, we use the expression Ij44|l . In the first case 
every node gives a contribution different from zero and equal to each other, so with the same observations 
made for the Kagome lattice we find that the longitudinal conductivity is 

2 

(Txx = ^ (Af = = 0). (54) 
When M > and (j) — we find that the conductivity is 



19 1 ^hj + hl + M2 



2he 8 (27r)2 hl + hl + - 
19 1 M 



d^k 



2he 8 (27r)2 ^ hj + h'i + {hj + + + e^) 

but, as in the case of the Kagome lattice, this integral vanishes; so in this case the longitudinal conductivity 
is zero. 

The last case we consider is for M = ±3\/3at2 sin(0) ^ 0. Here the bands touch only in three points 
and these give a contribution to the conductivity since now there is a gap where before there were three 
nodes. Around the three nodes the term c is zero, so formally the problem is the same as that of the case 
M = and = and only in these three points the longitudinal conductivity is different from zero. The 
result is equal to half of what was found in the case in which M and are zero, that is 

2 

(Txx = ^ {M = ±3\/3at2 sin(0) ^ 0). (55) 
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Figure 10: Hall conductivity for electrons on the Kagome lattice. 



We can conclude that the longitudinal conductivity can be rewritten as 



9 

e TT 



with = 0, 1, 2 for the honeycomb lattice and fx = 0, 2 for the Kagome lattice. 



4 Discussion and Conclusions 

We have considornd both Haldano's model for electrons in a staggered flux on the honeycomb lattice 
and the model by Ohgushi, Murakami and Nagaosa for electrons in the presence of a canted spin-1/2 
texture on the Kagome lattice. We have shown how similar these two models are in that the transverse 
Hall conductivity a^y is quantized as ±e^/h as a function of the tuning parameter (e.g. the magnetic flux 
per plaquette, (j)). 

Whilst for the Hall conductivity axy we have obtained the same results both for the Kagome lattice 
in the presence of a spin texture (as found by Oshgushi et al. (2000)) and for the honeycomb lattice with 
staggered magnetic field (as found by Haldane (1988)), we stress that we have used a different method of 
calculation based on implementing the band structure of each model in the Kubo formula. Furthermore, 
we have explicitely evaluated in this way, and for the first time, also the longitudinal conductivity a^x 
starting from the Kubo formula. This quantity is also quantized in the absence of symmetry-breaking, 
non-ideal features of the system, but not in terms of integer multiples of e^/h. 

For the Kagome lattice model, we find metallic behavior for a ferromagnetic state of localized spins 
perpendicular to the plane of the lattice. This state has a vanishing flux {(p = 0) in each plaquette of 
the Kagome lattice. Metalhc behavior exists also for a canted state where the spins are inside the plane 
(6* = 7r/2, in this case the local flux is 4) = tt). The longitudinal conductivity is for both cases axx = e^Tr/4h 
and the Hall conductivity vanishes. In Fig. 10 we show the schematic behaviour of the Hall conductivity 
for the model defined on the Kagome lattice and as a function of the parameter (j). This is to be compared 
with the richer phase diagram for the Haldane model on the honeycomb lattice, reported in Fig. 11 also 
as a function of 4>. The longitudinal conductivity axx as evaluated in this work is shown schematically in 
Fig. 12. 

Removing some of the nodes in the DOS by breaking symmetries (like for the case of a square lattice 
with next-nearest neighbor terms) alters the Hall conductivity substantially. Also the introduction of 
disorder (e.g. slow fluctuations of the locahzed spins, fluctuations around the perfect canted spin texture) 
may remove some of the nodes and yield non-universal features in the transport properties. There is also 
another interesting effect due to disorder in our two-dimensional lattices. The longitudinal conductivity 
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Uxx is usually based on diffusion of cfiarge carriers; however, the diffusion coefficient D is infinite in 
our model, since there is no scattering in the absence of imperfections. Nevertheless, the longitudinal 
conductivity, expressed through the Einstein relation 



is finite thanks to a vanishing density of states at the nodes. The cancellation of the divergent diffusion 
coefficient and the vanishing density of states is subtle. Since there is scattering by impurities in a realistic 
system, a finite diffusion coefficient is more natural. On the other hand, impurities create additional states 
near the nodes such that a non-vanishing density of states exists. This effect was studied in the case of 
2D Dirac fermions with random scatterers (Ziegler, 1997, 1998, Ziegler and Jug, 1997). In particular, it 
was found that random scattering broadens the metallic state (Ziegler and .Jug, 1997), and the maximal 
conductivity value is lowered by a factor 1/(1 + (7/27r), where g is the strength of the random fluctuations. 
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Appendix A: Linear Response and Kubo Formula 

From the Kubo formula we know that the conductivity tensor can be written as (Madelung, 1978) 

<T^u = ^\im r e(''^+«)*Tr([po,r^]e-'^*/'^>e'^*/'^)di (57) 

where po is the Fermi function. Using the Green functions defined as 

G±{E) = {H/n + E±ir])-'^ (58) 

we can use the substitution 

e^'"'/^ = ± lim r e^'^^'G^iE)— it < 0) 



(59) 



the conductivity can be rewritten as 



/oo T 771 

Tr[po,r^]G^{E)j,G-{E + io)— (60) 



The current operator is 

j, = -^[H,r,] (61) 



Using this expression iteratively, we find that the conductivity can be rewritten as a sum of three terms 

= - Tr{poG+{E)j^G+{E)jM-{E + io)}^ 

~hj_ Tr{poG+{E)j,G.{E + oj)j^G.{E + u;)}— (62) 

e f°° ( 
+ -j Tr{poG+{E)\j,,r^]G.iE + Lo)}: 



dE 
2^' 
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